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Abstract 

We consider the Casimir force betweeen two dielectric bodies described 
by the plasma model and between two infinitely thin plasma sheets. In both 
cases in addition to the photon modes surface plasmons are present in the 
spectrum of the electromagnetic field. We investigate the contribution of 
both types of modes to the Casimir force and confirm resp. find in both 
models large compensations between the plasmon modes themselves and 
between them and the photon modes especially at large distances. Our 
conclusion is that the separation of the vacuum energy into plasmon and 
photon contributions must be handled with care except for the case of 
small separations. 



1 Introduction 



The Casimir force between macroscopic bodies is at present one of the most inter- 
esting macroscopic manifestations of quantum effects. Applications range from 
attempts to obtain new constraints on hypothetical fifths forces in high precision 
experiments to impact on micromachinery and nanotechnology. The influence of 
the real structure of the interacting bodies is of increasing importance. For a 
better understanding of the Casimir forces efforts are undertaken to investigate 
the contributions from separate parts of the spectrum of the electromagnetic 
field. For instance, it is assumed that that surface plasmons (waves propagating 
along the surface 1 ) dominate the Casimir force at small separations while at large 
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separations the photon modes, i.e., the waves propagating perpendicular to the 
surfaces dominate. 

In the present paper we investigate the Casimir effect for thin plasma sheets 
and for dielectric bodies described by the plasma model and can confirm the 
above assumption only partly. We find large cancelations between the plasmon 
modes and and between the plasmon and the photon modes at all distances. For 
the first model this was observed already in [2|. 

The considered models are: 

1. Two dielectric bodies described by the plasma model 

Actually we assume two half spaces filled with a dieelectric medium and 
separated by a plane gap of width L. The dielectric function is described 
by the plasma model, 

e (-) = l-^§, (1) 
where u p is the plasma frequency. 

2. Thin plasma sheets 

Barton considered in [SJ E] a model of an infinitely thin shell filled with 
a charged fluid in order to describe a single layer hexagonal structure of 
carbon atoms aimed to model Cqq or carbon nano tubes. The interaction of 
such plasma sheet with the electromagnetic field results in a set of matching 
conditions for the latter. We consider two parallel sheets at separation L. 

In fact, these models result in very similar formulas. Both can be described by 
matching conditions and the Casimir force is given by the Lifshitz formula or 
a simple modification of it. Both models have surface plasmons in their TM 
polarization. The first model is well known and nearly all quantities considered 
below can be found in literature, see e.g., [2] or [3] and references therein. We 
consider this model here in parallel to the second one which is less known in order 
to demonstrate the similarity between the two. 

In general, the Casimir force can be represented as resulting from the ground 
state energy of the electromagnetic field, 

(n) 

where the sum runs over all quantum numbers, i.e., over the whole physical spec- 
trum. In the considered cases it consists of the photon modes which propagate 
in the whole space. They are labeled by a wave vector k e M. 3 and a discrete 
index for the two polarizations, TE and TM. In addition for the TM polarization 
there are the surface plasmons which are characterized by a discrete imaginary 
component of the wave vector in the perpendicular direction, k z = in. 
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In general, internal degrees of freedom of the considered media, density os- 
cillations in the second model for instance, give their own contributions to the 
vacuum energy. These contributions are not considered in the present paper. 

Of course, Eq as given by Eq. (J2J) contains well known divergencies. These do 
not influence the Casimir force and can be removed for instance by subtracting 
the energy for infinite separation. After that (J2J is finite and represented by a 
converging sum/integral over photons and plasmons, 



However, the integrand in -Ephoton is highly oscillating and very inconvenient to 
handle. Therefor it is natural to move the integration path to the imaginary axis 
where the integrand becomes exponentially decreasing and where it is easy to 
perform the integration numerically. This is used in the Lifshitz formula too and 
gives directly E total . 

Below we calculate the vacuum energy £7tatai m this way for both models. For 
the first model we repeat known results, for the second model the calculation is 
probably new. Then we calculate separately the plasmon contributions which is 
an easy exercise. From this the photon contribution is then -Ephoton = -Eaotai — 
-Epiasmon and we can compare the relative weight of the individual contributions. 

Technically our calculation is based on the formalism introduced in jH] (see 
also in the review 0, section 5.1.1) where the vacuum energy is calculated in the 
background of a potential depending on one coordinate only. In this formulation 
the handling of bound states is particularly transparent. 

We start in the next section with a description of the models and a collection 
of the necessary formulas. In the third section we describe the formulas we use 
for the calculation of the vacuum energy. In the fourth section the calculations 
are carried out and the results are presented. Section 5 deals with the special case 
of two attractive delta potentials and the last section contains the conclusions. 

Starting from here we put h = 1 and consider c as the speed of light relative 
to the vacuum. 

2 Description of the models 

Let us start from a scalar field $ in a geometry of parallel planes. After Fourier 
transform in time and in the directions parallel to the planes we note the wave 
equation, 



where c is the speed of light which may depend on the medium. This equation 
is to be supplemented by the corresponding matching conditions at z = ±L02. 
In the following we need the scattering solutions to equation (jlj) which are those 




photon • 



(3) 
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solutions which have for z — > ±00 the asymptotic behavior 

t e lfc2 , (5) 

where r and i are the reflection and transmissision coefficients. The frequency of 
these solutions follows from the dispersion relation uo 2 /c 2 = k 2 + k 2 . 

A surface plasmon is a bound state in the one dimensional problem and ap- 
pears on the positive imaginary axis at k = in with a wave function after division 
by t 



(6) 



and t has a pole in k = %k whereas r0t is finite. 
The matching conditions are as follows. 



1. Dielectric bodies 

The field strengths components D± and Eu must be continuous across the 
surface. The polarizations of the electromagnetic field can be separated into 
transverse electric (TE) and transverse magnetic (TM) ones. The matching 
conditions for the appearing scalar functions are 

$+ - = 0, $' + - = 0, (TE) 

e+$ + - e_$_ = 0, $' + - = 0, (TM) 1 ' 

where $± and e± are the limiting values from the right and from the left. 
The dispersion with e(u), Eq. (JTJ), is 



u = Jul + k 2 + k 2 . (8) 

The reflection and scattering coefficients are well known. We note them 
first for one surface assuming e 7^ 1 on the left side and e = 1 on the right 
side. Because we have plane waves outside the surface we can write down 
the wave function directly, 

$( ki . z) -f e lkz + re~ lkz for z < 0, 
with uj 2 = k 2 + q 2 . Then the coefficients are 



2efc _ fc-eg 

k+eq ' k+eq' 



* = r=fe£ (™) 



(10) 
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2. Plasma sheets 

For the plasma sheet the matching conditions demand E\\ to be continuous 
across the surface and for the normal component AE Z = -jV||.£7|| to hold, 
where Q is a parameter inverse proportional to the density of carriers in 
the plasma sheet. In this model the polarizations can be separated into TE 
and TM too and the matching conditions read then 



$+-$_ = 0, $' + - = 2Q $, (TE) 
$' _ = 0, $ + - $_ = -24 (TM) 



The matching condition for the TE polarization is just the same as for 
a delta function potential on a plane with strength Q. The potential is 
repulsive because Q > follows from the physics of the model. As known, 
a repulsive delta potential does not accommodate bound states and so there 
are no surface plasmons in this model. However, below we will investigate 
the unphysical case of Q < too where bound states are present. The 
matching condition for the TM polarization is sometimes associated with 
a 5' potential although this is not correct (see for example [7], Appendix). 
The combination of parameters —2-% appears from the physics of the model 
and leads to a bound state resp. to a surface plasmons. In this model the 
dispersion is simply 

w = Jk* { + k 2 (12) 
and the reflection and transmission coefficients are 



* ~~ Lu2+ikn> r ~ J+ikcu- (TM) 



(13) 



For two parallel surfaces at distance L the transmission coefficient written in 
terms of the coefficients from single surfaces is (a formula of this kind can be 
found e.g. in jS]) 

= tit r exp(i(q - k)L) 
1 + tLrtr r e 2i <* L ' 

where t is the complex conjugated to t and the indices indicate the left resp. 
right surface. This matters only for the first model where the coefficients (|Tnj) 
are written for the left surface. That for the right one follow by interchanging k 
and q and substituting e — > 1/e. In this way for the dielectric model 

4kg 



t = ^fr+ff (TE) 

Mt+§) e2tqL 



Aekq 
jk+eg)' 2 



, fc+eg ) 



(15) 
(TM) 
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holds and for the two plasma shells we have 

-fc 2 

f _ (»fc-n) 2 (TE) 



1_( n V e2rt£ ' 

t = SgM . (TM) 

V^ 2 +jfifc/ 

The surface plasmons are given by the poles of t on the imaginary axis at k — in. 
For the dielectric model we get the equation 

^ = -oe~* (17) 
k — eq 



with e = 1 — ^ 2+ fcf _ K 2 an d = y A;jj ~ k2 - F° r cr = +1 the corresponding wave 
function is symmetric and for a = —1 it is antisymmetric. The antisymmet- 



ric plasmon exists only for k\\ > LO p /yl + oo p L/2. We denote the solutions of 
equation (fTTj) by K a (k\\,u p , L). We introduce a similar notation for the frequency, 

uj(k lb k) = yjul + k 2 + k 2 . (18) 

A single surface with e given by (JIJ has one plasmon with 



^single 



uJ 2 



+ \ hf + (19) 



which can be obtained from K a (kn, cu p , L) in the limit L — > oo. 

For the plasma shell model we note the corresponding formulas. For the TE 
polarization the equation for k is 

1 + ^ = -ae- kL . (20) 

It has solutions for Q < only and the antisymmetric solution has the additional 
constraint > 1/L. For one plane the solution is Single = — O. For the TM 
polarization the equation is 



UJ 2 — KUJ 



ae~ KL (21) 



with u) = yk 2 — k 2 . The solutions exist for Q > 0. For one plane the solution is 

^single — \ ( K / \J^ 2 + ~~ ^J- Here we also introduce notations with an explicit 
indication of the arguments, K a (kn,Q, L) and 



0J{ 



(k lu k) = Jk 2 + k 2 . (22) 



Finally we note that for the plasmons the inequalities 

k_ < Single < K + < u p (dielectric plasma model), 

K- < /tangle < k+ < ku (model of two plasma sheets) (23) 

hold. 

3 The vacuum energy 

The vacuum energy for background potentials depending on one coordinate had 
been considered repeatedly. We follow here the derivation given in jH], or equiv- 
alently, in [3], section 5.1.1, but we change some notations, s i2 — > r and s u — > t, 
for instance. The basic formula is 

1 f dkn . . f°° dk „ d , t(k) } 

E °=2jj^{^ k «'*^ + l ^^ k w a iH)f' (24) 

where t(k) is the reflection coefficient defined in the scattering problem (JHJ and 
the sum runs over the surface plasmons. This formula is quite general and a 
dependence of u and of t{k) on k\\ is allowed. By means of = e 2lS it is 
equivalent to representations in terms of scattering phase shifts or the mode 
density. 

In this formula the ultraviolet divergences are still present. We remove them 
by subtracting the doubled contribution from one single surface. This is equiva- 
lent to subtract the limit of large separation L. In the formulas for the transmis- 
sion coefficients, Eqs. (fT5j) and ([16)1 . this results in removing the factors in the 
numerators and in the plasmon contributions we substitute 

Uj(k\\, iK a ) — ► L0(k\\, Mo) - to(k\\, ^single)- (25) 

Obviously, we subtract distance independent terms which do not contribute to 
the Casimir force. 

By means of Eq. (j21| the vacuum energy is given as the sum (jSJ) of plas- 
mon and photon contributions and it consists of sum/integrals over the physical 
spectrum. While the (subtracted) plasmon contribution is fast converging, the 
(subtracted) photon contribution is highly oscillating. As already mentioned, it 
is sufficiently inconvenient to calculate it in this representation and we move the 
integration over k towards the imaginary axis. In that fast converging 

integral appears which can be evaluated numerically without big effort. This 
procedure is standard and well known. For instance, the integration over the 
imaginary axis is used in the Lifshitz formula although this is not always easy to 
see because of different choices of the variables in the integrations. 

The only nontrivial moment which deserves attention is the role of the surface 
plasmons. In [H] it was shown that in the process of deformation the integration 
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contour crosses poles of the logarithmic derivative in (J24J). These are just the 
poles which correspond to the bound states resp. to the surface plasmons. The 
crossings give extra contributions which just cancel the sum in (j24jl . As a result 
one obtains the representation 

1 f dku f°° dk „ ,J 

E " s E >«* = 2]wr- i -"<*''• **> m ln **> (26) 

for the tote/ vacuum energy. The lower integration boundary for k is ko = 

\Jup + fcjj for the dielectric bodies with u(k\\,ik) given by (fTHj) and A; = fc|| 

for the plasma sheets with u(k\\,ik) given by ()22)) . This formula can be viewed 
also as a generalization of the Lifshitz formula to the case when surface plasmons 
are present. 

Surface plasmons exist in both models in the TM polarization only. However, 
in the TE polarization in the plasma shell model, which is equivalent to delta 
function potentials on the planes, surface plasmons exist for a 'wrong' sign of the 
charge density parameter Q, i.e., for Q < 0. Although in the considered model 
this case is not physical it is interesting to consider it anyway. We do this in 
section |SJ 



4 Calculation of the total vacuum energy and 
of the plasmon contributions 



As shown in the preceding section the total vacuum energy is given by Eq. 
For purposes of numerical evaluation it is meaningful to integrate by parts to get 
rid of the derivative. After that the formula for the total vacuum energy reads 

1 f dku f°° dk k / • 7 \ 

2 J (2tt) 2 J ho 7T u{k\\,ik) 

and this is the form which we use in the present section. 

Now we collect the formulas for the transmission coefficients on the imaginary 
axis taking the subtractions into account. For the dielectric bodies 

m))- 1 = ' e ~ 2qL > ^ TE ) 

(tiik))' 1 = 1 - (^^] e~ 2qL , (TM) (28) 



k + eq t 

holds with q = ^k' 1 — u'* and e = 1 — uj2+l S_ k 2 and for the plasma sheets we have 



s 



(tmr = i-i qz^hm ) e ~ 2kL - (TM) (29) 

Inserting these formulas into (|2*7|) the calculation of the vacuum energy is an 
easy numerical exercise. For L — > oo in all cases the ideal conductor limes 
E = -7r 2 /240 L -3 is reobtained. 

Note that Eq.([27|) contains a double integral. In all cases one integration can 
be carried out explicitly. In both models for the TE polarizations this is trivial 
because t(ik) does not depend on ku. For the TM polarizations after changing 
the variables from ku to u the integration over uo can be carried out explicitly. 
This integration can be done by standard computer algebra programs and the 
result is a quite lengthy expression which we do not display here. Within the 
dielectric model the corresponding integration was done for the force, i.e., for the 
derivative of Eq. (|2*7jl with respect to the distance, in the appendix in 

The contribution of the surface plasmons follows from Eq. (|24jl and (|25jl . It 
reads 

^plasmon = ~J ^ { yj + ^ ~ «£ ( k\ | , U p , L) 



yJul + q-Kl^{k\\,u p ,L)) (30) 



for the dielectric bodies with the constraint ku > co p / a/1 + u> p L/2 for a = — 1 
and 

gasmen = \ J ^ ^ { yj h*\ ~ K* ( k\ | , fi, L) 

- xJ^-^ih^L)) (31) 

for the plasma sheets. All these quantities are shown in Fig. ^for the dielectric 
bodies as a function of the plasma frequency uj p and in Fig. El for the plasma 
sheets as a function of the parameter Q for a fixed distance L = 1. In both cases 
in the left panel we observe a significant compensation between the two plasmon 
contributions. In the right panel we see a compensation between the plasmon 
and the photon contributions. 

The vacuum energy can be represented in the form 

£totai = -T^Jjf^vL), (32) 

where the factor in front is the ideal conductor contribution and the dimensionless 
function f(u p L) describes the deviation from the ideal conductor in the considered 
model. In case of the plasma sheet model it depends on QL. For large argument 
it has the limit 

lim f(x) = 1. (33) 
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Figure 1: The vacuum energy of the TM polarization in the dielectric plasma 
model and its several constituents, observe the large compensation in both panels 




Figure 2: The vacuum energy of the TM polarization in the plasma sheet model 
and its several constituents, observe the large compensation in both panels 




Figure 3: The contributions from the TE and TM polarizations to the vacuum 
energy in the dielectric model (left panel) and in the plasma sheet model (right 
panel) normalized to the perfect conductor limes 
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For both models and both polarizations this function is shown in FigEl 

Taking into account the dependence of the vacuum energy on the parameters 
as given by Eq. (j32|) we see from the figures Fig ^ and Fig|2] that the relative 
compensation between the plasmons and the photon contribution increases with 
distance L. This can be discussed in more detail. We consider the individual 
plasmon contributions as given by Eqs. ()30j) and (}3~Tj) and represent these formulas 
in the form 

1 [°° 

-Eplasmon, a = 7- / dk\\ k\\ UJ a (uJ p , L, k\\) (34) 

with fcfj = Up/ a/1 + uo p L/2 for a = —1 in the dielectric model and A;Jj = 
otherwise. In (EH) we introduced the notation 



cv a (u p , L, k\\) = ^Jul + &jj - K a (k\\,u p ,L) - Ju^ + fcfj - K single (k\\,u p , L). (35) 

The K a (k\\,ujp, L) are given by Eqs. (|T7|) and (J2T|) and Single ( k\ , u p , L) by Eq. (fT9|) . 
Now we make the substitution ku — * ku/L in the integral in ()34j). 

1 f°° 

-Epiasmon, <t = -r^r / dku ku Luj a (uj p , L, fen/ L) . (36) 

^ L Jk°L 

For dimensional reasons the expression Lu a (u p , L, ku/L) depends on the product 
ujpL only. This is a representation in parallel to ()32j) and the integral in ()36|) is 
a dimensionless function on u p L. In the next step we substitute ku —> k\\^Ju) p L 
and represent the plasmon contribution in the form 

ujpL 

^plasmon, a = T/~73 di^p^) (37) 

with 

g(u) p L) — y/(JpL / dk\\ k\\ Luj a {LOp,L,k\\Juj p /L). (38) 

So far this is a simple rewriting of the plasmon contribution. Now the statement 
is that 

• for the dielectric model the function g(x) has a finite limit for x — ► oo, 
namely 



lim g(:r) 



-1.2448 for a = 1 (symmetric plasmon), 
0.9652 for a = — 1 (antisymmetric plasmon). 

(39) 

The function g(x) is shown in Fig 0] for both plasmons. 



for the plasma sheets the function g(x) does not depend on its argument 
at all. The values are 

, , _ J —0.702427 for a — 1 (symmetric plasmon), , . 

1 0.639449 for a = — 1 (antisymmetric plasmon). 
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Figure 4: The function g(x) as defined in Eq.([38p in the dielectric plasma model, 
the lower curve is for the symmetric plasmon and the upper curve is for the 
antisymmetric plasmon 

This is the behavior of the plasmon contributions for large both or one of two, 
distance L and plasma frequency uj p (in the second model uj p must be substituted 
by VL). We see that the plasmon contributions grow by a factor ^Jui p L faster 
than the total vacuum energy. At once, by means of Eq.(j2J), this gives us the 
asymptotic behavior of the photon contribution and we see how the mentioned 
compensation is growing. 

It should be mentioned that in this way we got at once in the plasma sheet 
model the behavior of the plasmon contributions for small fl p and for small L. 

We note that the statements of Eqs. (fHHJ) and (|4f)|) are obtained numerically. 
No attempt had been made to get these results analytically. For the dielectric 
plasma model similar results are reported in [2]. 

In a similar manner the asymptotic behavior for small separations can be 
obtained numerically. In the dielectric model we find 



E TE ~ -0.002 a;;, 



E TM ~ -0.0039 

-^plasmon, cr=+l ~ —0.0305 — , 

^plasmon, .=+1 ~ 0.0267 (41) 

This confirms the known behavior that in this limit the TE contribution is small 
and the TM contribution is the sum of the two plasmon contributions. 
In the plasma sheet model we find 

Q 2 

E TE ~ -0.013 — , 

Lj 

E TM 0.00501 p5 
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For the plasmons, by means of Eq. (|40p. we know the behavior for all distances, 

-^plasmon, <r=+l 
-^plasmon, cr=— 1 

So in this model we observe the same pattern, the TE contribution is small and the 
TM contribution is at small separation the sum of the two plasmon contributions. 



0.05090 v^. 



-0.05589 -^r. (42) 
TP/- 1 



5 Two attractive delta potentials 

As mentioned above for the wrong sign of Q in the plasma sheet model the TE 
polarization has surface plasmons too. This is equivalent to a scalar problem with 
attractive delta potentials on two planes. An single attractive delta potential has 
always one bound state. For two planes there are two of them, provided L > Q 
holds otherwise there is only one. Further we note that, because the considered 
quantum field is massless, any bound state makes the potential instable, i.e., 
it causes particle creation. As a consequence, the vacuum energy acquires an 
imaginary part. In general, there is nothing special about that and we include 
this case here only because we are discussing the role of the surface plasmons. 

There are two ways to perform the calculations. The first is to make an ana- 
lytic continuation in Eq. (f2*U|) . When moving the binding energy K a of a plasmon 
which by means of (|23|) is initially below the integration part upwards so that k g 
becomes larger than k than one needs to deform the integration path in order it 
not to hit the pole. Than the integral remains finite but it acquires an imaginary 
part. 

The second way is to calculate the plasmon and photon contributions sepa- 
rately. Than it is obvious that the photon contribution is real and the plasmon 
contribution delivers just the imaginary part. Indeed, consider a plasmon contri- 
bution, 

-^plasmon = ~ J ^j) 2 ^ 7 

where k is a solution of Eq. (J20j) . As usual, we are concerned with a divergent 
quantity. Here the subtraction (j23|) of the contributions of single planes does not 
work because k does not depend on k\\. So we take a standard regular izat ion, say 

E* =1 [JtL A fkj^ C -V*IR (44) 

-^plasmon ^ J (^n) 2 V 11 ' 

This integral can easily be calculated explicitly, 

gasmen =^ + ^ + 0(5). (45) 
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After dropping the distance independent divergent part we just get 



7—1 / A r>\ 

^plasmon , i \*®) 

VZtx 

which is the contribution of one surface plasmon for two attractive delta poten- 
tials. 

In the photon contribution we proceed in the same way as before. We move 
the integration path for k and contract it on the imaginary axis. Now we have 
for k\\ < k a pole on the cut. It can be shown, that after the contraction of the 
path, for example using the well known formula 

1 1 

= Vp- +2vr<5(x), 

x — i\) x 

the imaginary part disappears and only the Value principal integral remains which 
is real. Actually this happens because the factor (k 2 + fcjj)2 has different signs on 
both sides of the cut. After that the order of integrations can be interchanged 
and the integration over k\\ can be carried out. Finally, integrating by parts, the 
photon contribution can be written in the form 

photon = h(2QL) (47) 



(4$ 



with 

45 f°° 
hh) = -r Vp/ dk k 2 log 

The function h(x) is shown in FigJSJ It has the property lim x -+± 00 h{x) = 1. For 
x > it represents the deviation of the vacuum energy of the TE polarization in 
the plasma sheet model from the ideal conductor limit and for x < it represents 
the deviation of the real part of the vacuum energy of two attractive delta poten- 
tials. We observe two regions with an attractive force and one with a repulsive 
force. 



6 Conclusions 

We calculated the vacuum energy in two models involving surface plasmons. As a 
result we are able to support the physically instructive picture that at small dis- 
tances the plasmon contributions dominate the vacuum energy. However, there 
is a compensation by one order of magnitude between the two plasmon contribu- 
tions (see Eq. (j41J) and (02)). At large separation we find that the compensation 
between the plasmon contributions grows and that in addition there is a large 
and growing compensation between the sum of the two plasmon contributions 
and the photon contribution. 
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Figure 5: In the left panel the function h(x), ([48)1 . which describes the vacuum 
energy relative to the ideal conductor limes for two delta potentials which are 
attractive for negative x and repulsive for positive x where they coincide with 
the TE polarization in the plasma sheet model, in the right panel the same 
function for small x 

In the dielectric model the compensations were found in the paper j^j where 
the importance of the plasmon contributions was emphasized. We would like to 
mention that the vacuum energy which is a relatively small quantity as compared 
to the plasmon contributions (even for for small separations there is a compen- 
sation between the individual plasmons, see Eq.fjUJ) can be directly calculated 
in the known representation where the momentum integration in the direction in 
perpendicular to the surfaces goes over the imaginary axis. In that formula large 
compensations are absent. A possible disadvantage, however, is that one needs 
to know e(u>) for imaginary uj. 

Also we considered the example of two planes carrying attractive delta poten- 
tials. Here the surface plasmons do not contribute to the real part of the vacuum 
energy. They deliver the imaginary part which is present because in that model 
particle creation is present. It is obvious that the imaginary part is present for 
large separations too. 

In general, the representations of the vacuum energy by an integral over the 
imaginary axis on the one hand side and as a sum of plasmon and photon contri- 
butions are connected by the dispersion relation the transmission coefficient t(k) 
obeys, see for example Eq.(27) in 6\. To some extend this is equivalent to the 
Kramers-Kronig relation for e(u). Both are exact relations and therefore a large 
compensation does not present a problem. However, if trying to make some kind 
of approximation in these relations the compensation may become a source of 
large errors. 
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